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We examine the validity of the results obtained with the singularity confinement integra- 
bility criterion in the case of discrete Painleve equations. The method used is based on 
the requirement of non-exponential growth of the homogeneous degree of the iterate of 
the mapping. We show that when we start from an integrable autonomous mapping and 
deautonomise it using singularity confinement the degrees of growth of the nonautonomous 
mapping and of the autonomous one are identical. Thus this low-growth based approach 
is compatible with the integrability of the results obtained through singularity confine- 
ment. The origin of the singularity confinement property and its necessary character for 
integrability are also analysed. 



The singularity confinement property has been proposed some years ago [1] as a discrete 
integrabihty criterion. The essence of the method is the observation that in integrable 
mappings a spontaneously appearing singularity does not propagate ad inGnitum under 
the action of the mapping but disappears ( "is confined" ) after some iteration steps. Thus 
singularity confinement appeared as a necessary condition for discrete integrabihty. How- 
ever the sufficiency of the criterion was not unambiguously established. The attitude (of 
the present authors at least) has always been that if the singularity confinement condi- 
tion were strong enough then it would suffice for integrablity, in perfect analogy with 
the Painleve-ARS [2] property for continuous systems. This sufficiency of the singularity 
confinement criterion was recently challenged by Hietarinta and Viallet [3] who produced 
explicit examples of mappings satisfying singularity confinement which are not integrable 
to the point of exhibiting chaotic behaviour. Their approach is based on the relation 
of discrete integrabihty and the complexity of the evolution introduced by Arnold and 
Veselov. According to Arnold [4] the complexity (in the case of mappings of the plane) 
is the number of intersection points of a fixed curve with the image of a second curve 
obtained under the mapping at hand. While the complexity grows exponentially with the 
iteration for generic mappings, it can be shown [5] to grow only polynomially for a large 
class of integrable mappings. As Veselov points out, "integrabihty has an essential corre- 
lation with the weak growth of certain characteristics" . Thus the authors of [3] proposed 
to directly test the degree of the successive iterates and introduced the notion of algebraic 
entropy. The method is appropriate for birational mappings. One starts by introducing 
homogeneous coordinates and studies the degree of the iterate. As Bellon and Viallet [6] 
remark, the growth of the degree is invariant under coordinate changes though the degree 
itself is not. A generic (non integrable) mapping leads to degrees that grow exponentially. 
The algebraic entropy is thus naturally defined a,s E = lim^^oo ^og{dn)/n where dn is the 
degree of the n-th iterate. Thus nonintegrable mappings have nonzero algebraic entropy. 
The conjecture in [3,6] is that integrabihty implies polynomial growth, leading to zero 
algebraic entropy. 

The main application of the singularity confinement approach was the derivation and 
study of discrete Painleve equations (d-P's) [7]. On the light of the results of Hietarinta 
and Viallet which have shown that the criterion used was not restrictive enough, one might 
be tempted to doubt the integrabihty of the mappings obtained (despite a considerable 
volume of integrability-confirming results). The aim of this paper is to show that these 
doubts are unjustified and to confirm the validity of the approach previously used, with 
the help of algebraic entropy techniques. 

Let us first recall what has always been our approach to the derivation of d-P's. We start 
from an autonomous system the integrabihty of which has been independently established. 



In the case of d-P's, this system is the QRT mapping [8]: 

f^'\xn) - (Xn+l + Xn-l)f^^\xn) + Xn+lXn-lf^^'Hxn) = (1) 

When the /(*)'s are quartic functions, satisfying specific constraints, the mapping (1) is 
integrable in terms of eUiptic functions. Since the eUiptic functions are the autonomous 
hmits of the Painleve transcendents, the mapping (1) is the appropriate starting point 
for the construction of the nonautonomous discrete systems which are the analogues of 
the Painleve equations. The procedure we used, often referred to as 'deautonomisation', 
consists in finding the dependence of the coefficients of the quartic polynomials appearing 
in (1) with respect to the independent variable n, which is compatible with the singularity 
confinement property. Namely, the n-dependence is obtained by asking that the singular- 
ities are indeed confined. One rule that has always been used, albeit often tacitly, is that 
confinement must be implemented "the soonest possible" . What this rule really means is 
that the singularity pattern of the deautonomised mapping must be the same as the one 
of the autonomous mapping. Our claim is that a deautonomisation with a different singu- 
larity pattern (for instance a 'later' confinement) would lead to a non-integrable system. 
The reason why this deautonomisation procedure can be justified is the following. Since 
the autonomous starting point is integrable, it is expected that the growth of the degree 
of the iterates is polynomial. Now it turns out that the application of the singularity 
confinement deautonomisation corresponds to the requirement that the nonautonomous 
mappings lead to the same factorizations and subsequent simplifications and have pre- 
cisely the same growth properties as the autonomous ones. These considerations will be 
made more transparent thanks to the examples we present in what follows. 

Let us start with a simple case. We consider the mapping: 

ClXji -\- , , 

Xn+1 + Xn-l = 2 ^^> 

where a and b are constants. In order to compute the degree of the iterates we introduce 
homogeneous coordinates by taking xo=p, xi=q/r, assuming that the degree oi p is zero, 
and compute the degree of homogeneity in q and r at every iteration. We could have of 
course introduced a different choice for xq but it turns out that the choice of a zero-degree 
xq considerably simplifies the calculations. We obtain thus the degrees: 0, 1, 2, 5, 8, 13, 
18, 25, 32, 41, . . . , . Clearly the degree growth is polynomial. We have d2m = 2m^ and 
d2m+i = 2m^ + 2m + 1. This is in perfect agreement with the fact that the mapping 
(2) is integrable (in terms of elliptic functions), being a member of the QRT family of 
integrable mappings. (A remark is necessary at this point. In order to obtain a closed- 
form expression for the degrees of the iterates, we start by computing a sufficient number 



of them. Once the expression of the degree has been heuristicaUy estabhshed we compute 
the next few ones and check that they agree with the analytical expression predicted). 
We now turn to the deautonomisation of the mapping. The singularity confinement result 
is that a and b must satisfy the conditions an+i — 2an + a^-i = 0, bn+i = &n-i, i-e. a 
is linear in n while 6 is a constant with an even/odd dependence. Assuming now that 
a and b are arbitrary functions of n we compute the degrees of the iterates of (2). We 
obtain successively 0, 1, 2, 5, 10, 21, 42, 85,.... The growth is now exponential, the 
degrees behaving like d2m-i = (2^"^ — l)/3 and d2m = 2(i2m-i7 a clear indication that the 
mapping is not integrable in general. Already at the fourth iteration the degrees differ in 
the autonomous and nonautonomous cases. Our approach consists in requiring that the 
degree in the nonautonomous case be identical to the one obtained in the autonomous 
one. If we implement the requirement that d^ be 8 instead of 10 we find two conditions 
an+i — 2an + an-i = 0, bn+i = bn-1, i.e. precisely the ones obtained through singularity 
confinement. Moreover, once these two conditions are satisfied, the subsequent degrees of 
the nonautonomous case coincide with that of the autonomous one. Thus this mapping, 
leading to polynomial growth, should be integrable, and, in fact, it is. As we have shown in 
[9], where we presented its Lax pair, equation (2) with a{n) = an + (3 and b constant (the 
even-odd dependence can be gauged out by a parity-dependent rescaling of the variable 
x) is a discrete form of the Painleve I equation. In the examples that follow, we shall 
show that in all cases the nonautonomous form of an integrable mapping obtained through 
singularity confinement leads to exactly the same degrees of the iterates as the autonomous 
one. 

Our second example is a multiplicative mapping: 

where one can put 6 = 1 through an appropriate gauge. In the autonomous case we 
obtain, starting with xq=p and xi=q/r, successively the degrees: 0, 1, 2, 3, 4, 7, 10, 
13, 16, 21, 26, . . ., i.e. again a quadratic growth. In fact, if n is of the form Am + /c, 
(A;=0,l,2,3) the degree is given by d^ = Aw? -f- (2m -|- l)k. The deautonomisation of (3) 
is straightforward. We compute the successive degrees and find: 0, 1, 2, 3, 4, 7, 11, ..., 
. At this stage we require that a factorization occur in order to bring the degree d^ from 
11 to 10. The condition for this is an+2ttn-2 = ^n i-^- '^ of ^^e form ae,oAg „ with an 
even-odd dependence which can be easily gauged away. This condition is sufficient in 
order to bring the degrees of the successive iterates down to the values obtained in the 
autonomous case. Quite expectedly the condition on a is precisely the one obtained by 
singularity confinement. The Lax pair of (3) can be easily obtained from our results in 



[10]. We find that if we introduce tlie matrices: Lr 
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we can obtain from tlie compatibility L„_|_iM„(/i/q') 



Mn{h)Ln tfie equation Xn+iXn-i = knkn+i{xn + l)/a^n7 wliere kn+i = qkn-i, wliicli is 
equivalent to (3) up to a gauge transformation. 

Tfie case of the mapping 

Xn+l + Xn-1 = an + -^ (4) 



has a more interesting deautonomisation. In this case we make a slightly different choice 
of homogeneous coordinates, which simplifies the results for the degrees of the iterates. 
We assume xq = p/r, xi=q/r and compute the degree of homogeneity in p, q and r. We 
find dn=l, 1, 2, 3, 5, 8, 11, 15, 20, 25, 31, 38, . . . , i.e. if n is of the form 3m we have 
dn='inn? — TO + 1, for n=3TO+l, dn=?>w? + m + 1, and for n=3m+2, dn='iTin? + 3m + 2. 
In the generic nonautonomous case the corresponding degrees are 1, 1, 2, 3, 5, 8, 13, . . . ,. 
The requirement that dQ=ll leads to the condition an+i = a^-i and 6n+2 — ^n+i — ^n-i + 
6n-2=0. Thus h is linear with a ternary symmetry while a is a constant (with an even/odd 
dependence which can be gauged away) . This fully nonautonomous form of (4) is a discrete 
form of Painleve IV studied in [11] and [12] where we have given its Lax pair. 

We now turn to what is known as the "standard" discrete Painleve equations [7] and 
compare the results of singularity confinement to those of the algebraic entropy approach. 
We start with d-Pi in the form: 

Xn+l+Xn+Xn-l=an + -^ (5) 

X"n 

The degrees of the iterates of the autonomous mapping are 0, 1, 2, 3, 6, 9, 12, 17, 22, . . . , 
i.e. a quadratic growth with (i3^+jt=3m^ + (2m+l)fc, for k =0,1,2 while those of the generic 
nonautonomous one are 0, 1, 2, 3, 6, 11, . . . , . Requiring two extra factorisations at that 
level (so as to bring d^ down to 9) we find the following conditions an+i = a^, so a must be 
a constant, and 6tj,+2 — &n+i — ^n + &n-i=0, i.e. 6 is of the form 6^ = Q;n + /3 + 7(—l)"' which 
are exactly the result of singularity confinement. Implementing these conditions we find 
that the autonomous and nonautonomous mappings have the same (polynomial) growth 
[6] . Both are integrable, the Lax pair of the nonautonomous one, namely d-Pi having been 
given in [10,13,14]. 



For the discrete Pn equation we have 

Xn+1 + Xn-1 = ^——. (6) 

X^ i 

The degrees of the iterates in the autonomous case are dn=0, 1, 2, 4, 6, 9, 12, 16, 20, 
. . . , (i.e. (i2m-i="^^5 d2m='n^'^ + ^) while in the generic nonautonomous case we find the 
first discrepancy for (i4 which is now 8. To bring it down to 6 we find two conditions, 
ttn+i — 2an + an-i=0 and 6n+i = ^n-i- This means that a is hnear in n and h is an 
even/odd constant, as predicted by singularity confinement. Once we implement these 
constraints, the degrees of the nonautonomous and autonomous cases coincide. The Lax 
pair of equation (6) in the nonautonomous form, i.e. d-Pn, has been presented in [10,15,16]. 

The Q-Piii equation was obtained from the deautonomisation of the mapping: 

Xn+lXn—1 = '77~2_ \r\~Z. / V ^ ' 

In the autonomous case we obtain the degrees (in=0, 1, 2, 5, 8, 13, 18, . . . , just like for 
equation (2), while in the generic nonautonomous case we have 0, 1, 2, 5, 12,..., . For (i4 
to be 8 instead of 12, one needs four factors to cancel out. The conditions are c^+i = c„_i 
and a„_|_i6„_i = an-ihn+i = anbn- Thus c is a constant up to an even/odd dependence, 
while a and b are proportional to A"^ for some A, with an extra even/odd dependence, just 
as predicted by singularity confinement in [7]. The Lax pair for q-Pm has been presented 
in [10,17]. 

For the remaining three discrete Painleve equations the Lax pairs are not known yet. It 
is thus important to have one more check of their integrability provided by the algebraic 
entropy approach. We start with d-Piv in the form: 

{xn+i + Xn)(:r^-i + Xr,) = ^ - '\^^_ ^ ' (8) 

\Xn \ ZjiJ C 

where a, b and c are constants. If z^ is constant we obtain for the degrees of the successive 
iterates (in=0, 1, 3, 6, 11, 17, 24, . . . , . The general expression of the growth is dn=Qm'^ 
if n = 3?B, dn=Q'w? + Am + 1 if n = 3?ti + 1 and dn=Qni? + 8m + 3 if n = 3m + 2. This 
polynomial (quadratic) growth is expected since in the autonomous case this equation is 
integrable, its solution being given in terms of elliptic functions. For a generic Zn we obtain 
the sequence (in=0, 1, 3, 6, 13, . . . , . The condition for the extra factorizations to occur in 
the last case, bringing down the degree (i4 to 11, is for z to be linear in n. We can check 
that the subsequent degrees coincide with those of the autonomous case. 

For the Q-Py we start from: 

/ ^ _ 1 V'T 'T _ 1 "l — l-^n ~r (i'Xn + -l-Jl,^:^^ ~r '^•^n + -l-j /„n 

[Xn-'rIXn ^)\Xn — lXn J-J /^ x /^ , x \d ) 

\ J- Z'YiCXjijyi. ZjiCiXjij 



where a, 6, c and d are constants. If moreover z is also a constant, we obtain exactly 
the same sequence of degrees dn=0, 1, 3, 6, 11, 17, 24, . . . , as in the d-Piv case. Again, 
this polynomial (quadratic) growth is expected since this mapping is also integrable in 
terms of elliptic functions. For the generic nonautonomous case we again find the sequence 
dn=0, 1, 3, 6, 13, . . . , . Once more we require a factorization bringing down (i4 to 11. It 
turns out that this entails a z which is exponential in n, which then generates the same 
sequence of degrees as the autonomous case. In both the d-Piv and q-Py cases we find 
the n-dependence already obtained through singularity confinement. Since this results to 
a vanishing algebraic entropy we expect both equations to be integrable. 

The final system we shall study is the one related to the discrete Pyi equation: 

where a, 6, c and d are constants. In fact the generic symmetric QRT mapping can be 
brought to the autonomous {z^ constant) form of equation (10) through the appropriate 
homographic transformation. In the autonomous case, we obtain the degree sequence 
dn=0, 1, 4, 9, 16, 25, . . . , i.e. dn=n^ ■ Since mapping (10) is rather complicated we cannot 
investigate its full freedom. Still we were able to perform two interesting calculations. 
First, assume that in the rhs instead of the function Zn a different function C^ appears. In 
this case the degrees grow like 0, 1, 5, ... , and the condition to have (i2=4 instead of 5 
is Zn-s^iZn-iz^ = (^. Assumiug this is true, we compute the degree d^ of the next iterate 
and find (is =13 instead of 9. To bring down ds to the value 9 we need z'^ = Cn^ which up 
to a redefinition of a and b means Zn = Cn- This implies Zn+iZn-i = z'^, and Zn is thus an 
exponential function of n, 2^=A" (which is in agreement with the results of [18]). Then a 
quartic factor drops out and ds is just 9. One can then check that the next degree is 16, 
just as in the autonomous case. Thus the q-Pyi equation leads to the same growth as the 
generic symmetric QRT mapping and is thus expected to be integrable. As a matter of 
fact we were able to show that the generic asymmetric QRT mapping leads to the same 
growth dn=n'^ as the symmetric one. This is not surprising, given the integrability of this 
mapping. What is interesting is that the growth of the generic symmetric and asymmetric 
QRT mappings are the same. Thus dn='n? is the maximal growth one can obtain for the 
QRT mapping in the homogeneous variables we are using. As a matter of fact we have 
also checked that the asymmetric nonautonomous Q-Pyi equation, introduced in [18] led 
to exactly the same degree growth dn='n? . 

Let us summarize our findings. In this paper, we have compared the method of singular- 
ity confinement and the approach based on the study of algebraic entropy when applied 
to the deautonomisation of integrable mappings. We have shown that in every case the 



confinement condition which ensured that the singularity pattern of the autonomous and 
non-autonomous cases are identical was precisely the one necessary in order to bring the 
growth down to the one obtained in the autonomous case. This validates the deautonomi- 
sation results obtained through singularity confinement at least in the domain of d-P's. 
This suggests also a strategy for the study of integrable mappings. We believe that in 
the light of the present results, when one starts from an integrable autonomous mapping, 
the deautonomisation can be performed solely with the help of singularity confinement, a 
procedure considerably simpler than the calculation of the algebraic entropy. 

Our present investigation also sheds light on the singularity confinement, and its necessary 
character as discrete integrability criterion. Let us go back to the example of mapping (2) 
with 6=1. We start with xq = p, xi = q/r. Iterating further we find 

r^ + aqr — pq^ qP4 (r^ + aqr — pq^)PQ -P4-P9 

where the PkS are homogeneous polynomials in q, r of degree k. (Remember that p is of 
zero homogeneous degree in our convention). The pattern now becomes clear. Whenever a 
new polynomial appears in the numerator of Xn its square will appear in the denominator 
of a:^_|_i and it will appear one last time as a factor of the numerator of x^_|_2 7 after 
which it disappears due to factorisations. The singularities we are working with in the 
singularity confinement approach correspond to the zeros of any of these polynomials, 
which explains the pattern {0, oo^,0}. The singularity confinement is intimately related 
to this factorisation which plays a crucial role in the algebraic entropy approach. Let us 
suppose now that a is a generic function of n. In this case we get the sequence: 

r'^ + aiqr — pq^ qQ4 {r"^ + aiqr — pq^)Qi 

qz ryr'^ + aiqr — pq^y qQ 
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qQ^Qi2 



r{r'^ + aiqr — pq'^)Q^ 

where the Qfc's are also homogeneous polynomials in g, r of degree k. Now the simplifica- 
tions that do occur are insufficient to curb the asymptotic growth. As a matter of fact, if 
we follow a particular factor we can check that it keeps appearing either in the numerator 
or the denominator (where its degree is alternatively 1 and 2). This corresponds to the 
unconfined singularity pattern {0, 00^, 0, 00, 0, 00^, 0, 00, . . .}. Once more, the confinement 
condition a^+i — 2a + a^-i = is the condition for q to divide exactly Q?, for both q and 
r^ -|- aiqr — pq^ to divide exactly Q12, etc.. Our analysis clearly shows why singularity 
confinement is necessary for integrability while not being sufficient in general. Still in the 
case of integrable deautonomisation it does lead to the correct answer, which explains its 
success in the derivation of the discrete Painleve equations. 
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